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Abstract: Inflation generically produces primordial gravitational waves with a red spec- 
tral tilt. In this paper we calculate the backreaction produced by these gravitational waves 
on the expansion of the universe. We find that in radiation domination the backreaction 
acts as a relativistic fluid, while in matter domination a small dark energy emerges with 
an equation of state WeS ~ —8/9. 
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1. Introduction 

The discovery of an apparent universal acceleration in the late 1990s [1, 2] lead to a mul- 
titude of theoretical explanations. Amongst the most appealing of these, given that they 
involve no new physics, are those based on the inhomogeneous structure of spacetime. 
Spatial inhomogeneities can influence cosmological models in two main ways: they alter 
the null geodesies and so affect photon propagation and the luminosity distance dL{z), 
and due to the nonlinearity of the Einstein tensor they cause the evolution of the average 
cosmology to deviate from the Robertson- Walker behaviour. In the language of [3], these 
are dubbed "weak cosmological backreaction" and "strong cosmological backreaction" re- 
spectively, although it might be noted that the term "backreaction" typically refers to the 
impact of perturbations on their background. The weak backreaction - deviations of the 
null geodesies from Robertson- Walker predictions - is generally studied with reference to 
LeMaitre-Tolman-Bondi (LTB) models (see for example [4, 5, 6, 7, 8, 9, 10, 11] and their 
references for a small sample of the recent literature) and we will not consider it further 
here. The so-called "strong backreaction" can be formulated in general coordinate systems 
[12, 13] and has a long history [14, 15, 16, 17, 18, 19, 20]. The discovery of the universal 
acceleration lead to an increase in interest in the field [21, 22, 23, 24, 25, 26]. Current re- 
search focuses on extending the formalism (see for example [12, 13, 27]) and on quantifying 
the impact in particular systems [28, 29, 30]. There has been much recent work in the field, 
and we point the interested reader to Buchert's 2007 review [31] and the references found 
therein and in [28] for further references. 

In this paper we focus on the backreaction arising in perturbed Friedmann-LeMaitre- 
Robertson- Walker (FLRW) universes, which serve as good models for the real universe 
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across a wide range of redshifts. The perturbations to this model can be separated into 
scalar, vector and tensor components based on how they transform on spatial three-surfaces. 
Of the three classes, the scalars have been widely considered since the earliest studies [14, 
15, 16] of cosmological averaging, since scalars are the dominant perturbation and develop 
strong inhomogeneities under structure formation. Marochnik [15, 16] also considered the 
vector and tensor modes, which have not received much attention since. Tensor modes 
decay as radiation and one would naively expect the backreaction they induce at the 
present day to be many orders of magnitude beneath that from the scalars - conversely, it is 
conceivable that, deep in radiation domination the tensor backreaction could be equivalent 
to or greater than that from scalars. The consideration of tensor modes is complicated 
by our lack of a simple covariant averaging procedure - the so-called Buchert approach, 
which employs a 3+1 spatial slicing and averages only scalars, cannot meaningfully be 
applied to tensor quantities. In principle, the covariant averaging approach of Zalaletdinov 
[32] can be applied to tensor modes, although attention with this approach has focused 
on solutions exhibiting a high symmetry [33, 34] or on scalar modes (e.g. [35]). Recent 
suggested averaging schemes by Behrend [36], Gromes [37, 38], Korzynski [27] and Coley 
[39] should also in principle be applicable to tensor modes. In this paper we will employ the 
standard 3+1 scalar averaging approach, focusing on the scalar quantities characterising 
the gravitational waves. 

Scalar backreaction produces multiple correction terms. We employ the terminology 
and approach of [40] (hereafter Paper I), [28] (hereafter Paper II) and [13] (Paper III). 
The corrections induced by the inhomogeneities are separated into six components: 7?.x), a 
correction to the Friedmann equation that arises from the average spatial curvature; the 
kinematical backreaction Qd, which appears in both the Friedmann and Raychaudhuri 
equations; the dynamical backreaction Vx>, which appears in the Raychaudhuri equation; 
and Tv, Sx> and Cj) which correct the fluid and cosmological constant properties. Paper II 
considered the scalar backreaction for a range of redshifts z G (0, 10^), finding that Tlx) is 
at least an order of magnitude greater than Vv-, while Qv is insignificant. This holds both 
in matter and in radiation domination. The effective equation of state of the backreaction 
evolves from w^q ~ 1/3 in radiation domination towards w^q ~ 1/60 at the present day^. 
In radiation domination, despite the effective equation of state, the effective energy density 
of the backreaction is approximately constant, while in matter domination it grows with 
the scale factor a. For a few recent studies focusing on scalar modes see [41, 40, 42, 28, 30]; 
the results are broadly consistent with one another when evaluated on the same scales, 
although the approaches are significantly different. 

In the absence of active sources such as primordial magnetic fields or cosmic defects, 
there are two main mechanisms by which cosmological gravitational waves are produced: at 
linear order, directly from inflation, and at second-order, sourced by the linear scalar modes. 
While one might naively expect the latter to be subdominant, they can instead at some 
periods and on some scales be of equivalent magnitude to the linear modes [43, 44, 45]. This 

^It should be noted that this figure is evaluated in Paper II assuming that averages of second-order 

perturbations of the form {4>{2)) vanish. This will not in general be true; see for example [30]. In Paper III 
it was argued that this remains approximately true for systems where {4\i)) ~ {'t>(2))- 
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could be particularly significant around the epoch of matter /radiation equality [44]. The 
analysis in an averaging context of such modes will be non-trivial due to the complicated 
form of their time-dependant power spectrum. In contrast, the statistics of the linear modes 
are imprinted at inflation and do not evolve. Furthermore, the equations of motion allow 
for analytic solutions, perhaps leading to closed forms for the backreaction they induce. 
We choose, therefore, to focus here only on the linear modes, giving us a good impression 
of the type of impacts gravitational radiation can have, while being amenable to analytic 
study. Modifications from the second-order modes would then be expected to be of at most 
the same order of magnitude, and we leave them to future study. 

As compared to the scalar case, one would perhaps expect the curvature correction TZt> 
induced by inflationary gravitational waves to remain dominant. However, there is no rea- 
son to believe that the tensor backreaction terms Qxi and Vv should remain subdominant 
or evolve in the same manner as those from the scalars; unlike the scalar perturbations, 
tensor modes are unaffected by structure formation. However, we would expect significant 
differences between the nature of the backreaction that arises in radiation domination and 
that which arises in matter domination; in radiation domination a first-order tensor mode 
evolves with conformal time r/ as ~ jo{kr]), while in matter domination it does so as 
h ~ hiji{kri) / (kr]) + h2yi{kri) / (krj) . This rapid decay in matter domination will induce 
differences between the evolutions of the correction terms TZd, Qt> and Vd, leading perhaps 
to distinctive differences between the effective energy density and pressure in the very early 
and in the recent universe. Since the tensor modes are not affected by structure formation, 
unlike the scalars the effective energy density of the backreaction in matter domination 
should not be expected to grow rapidly; while this suggests that it should be small in the 
present epoch, it conversely suggests there could be a period in the past when the tensor 
backreaction dominates the scalar. The two main aims of this paper are to consider these 
problems: firstly, to determine the impact of the tensor modes deep in radiation domination 
and whether it is ever plausible that it is comparable to the scalar modes; and secondly, 
to investigate the behaviour of the tensor backreaction when the tensor modes decay in 
matter domination, and the nature of it today. 

In §2 we set up the averaging formalism which we will employ, and in §3 we apply it 
to an FLRW universe containing scalar and tensor perturbations. §4 presents our findings 
and we end with a brief discussion in §5. 

2. Formalism 

Let spacetime be foliated on large scales with a family of 3-surfaces, choosing a coordinate 
system with a vanishing shift vector and the line element 

ds'^ = -a^dt^ + hijdx'dxK (2.1) 

The surfaces have the normal 4-vector 

n'^ = (l/a,0), n^ = (-a,0) (2.2) 
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with the projection tensor onto the surface 

hfj,u = g/xu + n^rit,. (2.3) 

The energy density and isotropic stress on this surface are 

Q = ni^n^T^,, S = h'^Sij = h'^h^h^T^, = hi^'T^, (2.4) 

and are related to the rest-frame energy density p and isotropic pressure p of a fluid with 
4-velocity by 

Q = p+ in'^n" - u^u^) r^^, S = 2p+ {n^'n^ - u^^u^) T^^. (2.5) 

Following for example Paper III we define a volume, average and an average volume ex- 
pansion in a domain V lying on the 3-surface by 

V-n= I Vhd^^, {A) = ^ [ AVhd'^, s^ = ^ (2.6) 

where h is the determinant of the spatial 3-mctric. Then volume- averaged Fricdmann and 
Raychaudhuri equations can be derived from the Hamilton constraint and evolution of the 
extrinsic curvature respectively, 

(S) ' = ^ («^P> + \ {o?) ^-\{nv + Qv- 6Tv) , (2.7) 

^ = (P + 3p)> + \ {a') A + Uvv + Qv- 3^2,) . (2.8) 

Here the curvature correction, kinematical backreaction, dynamical backreaction, and tilt 
correction are 

Uv = {a^n) , (2.9) 
Qv = {a' {k' - KiKl) ) - ^ {aKf , (2.10) 
Vv = {aV'Via) - {aK) , (2.11) 
-^i? = ^ <«' (n'^n^ - u^'un T^,) . (2.12) 

Here Kij is the extrinsic curvature and H the Ricci scalar on the 3-surface. These expres- 
sions are valid for any system for which the 3-1-1 split holds. 

3. Tensor Perturbations 

In Paper III the general backreaction expressions in a spatially flat Robertson- Walker 
universe perturbed to second-order with line element 

ds"^ = a^{ri) (- (1 + 20) drf + 2Badr]dx'' + {6ij + 2Cij) dx'dx^) (3.1) 
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were presented. We wish to consider a universe containing scalar and tensor perturbations, 
neglecting the vector modes which, at linear order, decay rapidly. It is important to 
note that, while the tensor modes at linear order are intrinsically gauge-invariant, gauge 
issues remain a problem when considering their backreaction for two chief reasons. Firstly, 
the spatial volume is dependant upon ^fh and, secondly, the tensor modes will couple 
to gradients of the scalar modes. We must therefore specify a gauge in which to work. 
Although a strong argument can be made for working in flat gauge [13], we work for 
simplicity in Poisson gauge, which takes the line element 

ds^ = a^ir]) (- (1 + 2(1)) drf + ((1 - 2il))5ij + hjj) dx'dx^) (3.2) 

where the tensor component of the 3-metric is gauge-invariant and transverse-traceless, 

S'^hjj = d'hfj = 0. (3.3) 

The same caveats highlighted in Papers II and III concerning the averaging of a perturba- 
tion theory also apply here. 

The curvature correction and the kinematical and dynamical backreactions are given 
by equations (40, 42, 43) of Paper III specialised to the case of a vanishing shift. Taking 

Cij = -i^dij + (3.4) 

and separating the corrections to the background into their scalar parts and tensor contri- 
butions, the curvature correction becomes 

TZv = + (^hp'^dkhfj + Ahf^d^d^^ + I [d%') {d,hj,) - \ (d^h^) {d,h[,)^ . (3.5) 

The scalar contribution TZ^^ is given in equation (82) of Paper III and does not concern 
us here. The kinematical and dynamical backreactions can likewise be written 

Qv = Q^^^-\{f^M)^ rv = vi,'^-(^h[,(d^d^<P + lhf^y (3.6) 

Their scalar components are given in equations (80) and (81) respectively of Paper III. For 
a perfect-fluid source the tilt correction J^d is unchanged by the presence of gravitational 
waves, as are the subsequent modifications Tij, 5x) and Cd which additionally correct for 
the non-unity lapse. The errors introduced by the coupling of the gravitational waves and 
the anisotropic stresses of the relativistic fluids is expected to be subdominant. The tensor 
corrections are intrinsically quadratic in that the leading-order backreaction arises from 
first-order squared terms and can thus be calculated good to second order in perturbation 
theory using linear perturbations and without complications from the metric determinant 
Vh. From (2.7, 2.8), the effective energy density and pressure of the backreaction can be 
identified as 

8^^„2- _ SttG 2-(S) ^(o^^^^V^^A f3 7l 

8^'^„2- _ SttG 2-(5) 1 /^(T) ^(T) .„(T)\ , . 
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In the remainder of this paper we wiU employ the normahsed variables 

= ^a^pi? = (3.10) 

eff 3^2 ^eff 18(d/a)2 ^ ^ 

and 

(T) ^ ^ (AR\ ^ 1 /^(T) , ^(T)^ 



eff 



n(I)' U J -3|(d/a) |r^ J ^'-^'^ 



^eff 

to characterise the contribution of the tensor modes to the background. These correspond 
to the effective energy density arising from tensor modes, the corresponding effective pres- 

(T) (T) 

sure, the effective equation of state, and the effective gravitational energy p^g + 3p^g con- 
tributing to the Raychaudhuri equation, all normalised to the background expansion. A 

(T) (T) 

positive ^^gfj increases the effective Hubble rate. An effective equation of state w^f^ < — 1/3 
is equivalent to a positive {AR/ R)^'^\ and in this case the tensorial backreaction will act 

(T) 

to accelerate the volume expansion. For an effective cosmological constant, w^^g' = —1 and 
(Ai?/i?)(-^) = 2r2gg\ while an effective dark energy with the equation of state WeS = —8/9 
would produce {AR/R)^'^^ = (5/3)^^^^. 

4. Results 

4.1 General Considerations 

The expressions (3.5, 3.6) are in general spatial integrations across real variables. The most 
frequent approach to these equations employs ergodic averaging (see for example [29, 30, 13] 
for recent discussions), working with domains large enough that a volume average tends 
towards the ensemble average. This approach has the advantage of being well-studied and 
amenable to both analytic and numerical work as the perturbations reduce to variables of 
the wavenumber only. 

Using the Fourier convention 

^(x) = J A(k)e-*'^-^, ^(k) = J A(x)e'»'-d3x (4.1) 

the average of a quantity quadratic in linear perturbations can be rewritten as 

<A2(x)) =a^l Wj,{^)A\^)d'^ = a^ll T^^(k)A(k')A*(k' - k)^^ (4.2) 

where W{x.) is a window function specifying the domain V, and 

Taking an ensemble average of (4.2) and employing the definition of the power spectrum 

^(k)^*(k') = ^\A{k)fv{k)6{k-k') (4.3) 
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we find that 
Here 

^ _ a3 J ^(x)d3x _ / W{^)^^d?^ ^ ^ 

Vv JW{x)^d^x ^ ^ ' ' 

is the ratio between the domain volume evaluated on the background and evaluated on the 
perturbed surface. The results of [40, 28, 42, 30] suggest that the averages of the scalar 
perturbations remain relatively small, implying that to a good approximation we can take 
i? ~ 1. It is now straightforward to show that terms containing contractions of a tensor 
mode with a partial derivative, such as (hj-d'' ijj^ , vanish upon ensemble averaging, since 
the tensor modes arc transverse-traceless. It is worth emphasising that, since the tensor 
modes are intrinsically gauge- invariant, the vanishing couplings with scalar quantities with 
i? 1 now implies that our results will themselves be automatically gauge-invariant. 
The tensor perturbation can be expanded into its two polarisation states by 

hl=h+Al. + h^A^.. (4.6) 

A^-'^ form an orthonormal basis and explicit forms could be found, for example, by align- 
ing the coordinate system with the Fourier mode (sec for example [46, 47]). The two 
polarisation states are assumed to possess similar statistics, that is that 



h+{k)h%{k') hy,{k)h*^{k') « h{k)h*ik) (4.7) 

We assume the primordial gravitational waves to have been produced by some inflationary 
process, with a power spectrum 

At is the amplitude of the primordial tensor perturbations, evaluated at a pivot wavenum- 
ber k^,, which is typically chosen to be /c* ^ 0.05Mpc~^. In the simplest inflationary models. 
At = rAs, where As is the amplitude of primordial scalar perturbations. The tensor spec- 
tral index is ut, which in simple models is approximately ut ~ — r/8. The gravitational 
wave spectrum is therefore slightly red and we will generally employ the inverse spectral 
index, m = —tit, which is a small, positive number. When a specific model is required we 
take r = 1/20 and m = —ut = 1/160. 

The tensor corrections to the Friedmann and Raychaudhuri equations (3.5,3.6) can 
now be written as 

= - (2^^ / ^rmmimf = 2,3'- [ kh'^, (4.10) 

irr\ IF ■ dk r°° ■ dk 
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with 

At linear order, the kincmatical backrcaction is related to the energy carried by the gravi- 
tational waves, while the dynamical backreaction is related to their energy flux. 
Gravitational waves evolve as 

h + 2-h + k'^h = ^ (4.13) 

where tt contains the fluid anisotropic stresses. In the following two sections we solve this 
equation and evaluate the backreaction terms (4.9-4.11) in radiation and matter domination 
respectively. 

4.2 Radiation Domination 

In radiation domination the photons and baryons are tightly coupled to one another and 
their anisotropic stresses driven towards zero. Neglecting the anisotropic stress of the 
neutrinos, we can set this term to zero. The gravitational wave in radiation domination is 
then 

h = (4.14) 

kr] 

where we have removed a decaying mode that otherwise diverges as ry — >■ 0. Using the 
coordinate x = kr], the curvature correction and backreactions (4.9-4.11) become 



?)^^fsin^.^, (4.15) 
,(T) ^W_r /sin2x _ sin^xX _dx_ 



(T) (3 [°° ( 2 sin 2a:; ^ sin^x^ dx 



Q}"' = -I— I ( cos^a;- ^^^^ + ^^^^ ) (4.17) 

For m > we can find closed-form exact solutions. The backreaction terms in radiation 
domination are 

(T) ^ /3 v/^ r(gf^) 

V (l + m)(2 + m)r(m)sin(7rm/2) (l + m)(2 + m) ^' 



(T) _ P 2-m TT _ 2 - ™ ^(T) 

V ~ ^2-m22-m(2 + 77i)r(l-Fm)sin(7rm/2) ~ 2-Fm ^' 



In the limit m <C 1 expected in a realistic model of inflation, these reduce to 

^^^^^2^^' ^ir^^-2m<\ Qg^)«(l-m)<^ (4.19) 

In Papers I and II it was shown that the dynamical backreaction arising from scalar per- 
turbations is an order of magnitude smaller than the scalar curvature, with the kinematical 
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backreaction negligibly small. In contrast, the kinematic backreaction produced by infla- 
tionary gravitational waves in radiation domination is approximately equal to the scalar 
curvature while the dynamic backreaction is suppressed by one power of the spectral index. 
The effective energy density (3.9) in radiation domination is given by 

= -l^rj'^J— -J-Pv"" (4.20) 

6^' 2 + m r(l + m) sin(7rm/2) 6m^ ' ^ ' 



and the effective normalised pressure (3.10) is 

t 

18 r(3 + m) sin(7rTO/2) ' 3 



Since ,0 < 0, both the effective energy density and pressure are positive definite. To linear 
order in m, the effective equation of state (3.11a) is 



« ^ (1 - 5m) , (4.22) 



demonstrating that in some respects the tensor backreaction fluid acts as an effective 
radiation. Finally, the normalised impact on the Raychaudhuri equation (3.11b) is 



Ai?\ _ 2"' (2 - 3;-/-/ - n?) „ tt _ 2 - 3/» - 



ni 



2 



«„™ - = '" o(^) (A 

Rj 12 ^ ' r(3 + m)sin(7rm/2) 2(1 + m) ' ^ ' 

As P < 0, (Ai?/i?)(^) < 0. The tensor backreaction in radiation domination then acts to 
decelerate the volume expansion. It is also worth noting that the effective fluid properties 

(T) 

grow only slowly with time, 0,^,^' (x rj"^ oc a"^ and so are almost constant. 
4.3 Matter Domination 

The situation in matter domination is more complicated. Continuing to assume a vanishing 
anisotropic stress is a good approximation while the photons remain tightly-coupled. Dur- 
ing and after recombination, the photon energy density is low enough that, even though the 
photon anisotropic stress is non-vanishing, the approximation remains at least reasonable. 
The solution to (4.13) in a universe dominated by dust is given by 

fa/cos(A:;/) - siu(A7/) /,:7/sin(/,:;/) +cos(fcr/) 
h = ^k + Bk . (4.24) 

The normalising functions and are found by matching the mode and its first deriva- 
tive at equality, and in terms oi x = kr] = ktrje are given by 

A, = -2 + cos^ (f ) 4 sin (I) , 5, = I + i sin (|) + 2^ sin^ (f ) . (4.25) 

Here t is the conformal time in units of the conformal time at matter /radiation equality, 
t = 1 corresponding to the epoch of equality and t 110 to the present epoch.^ The 



■^In deriving this value we have assumed flm ~ 1 and Ho = 70.1kms ^Mpc ^. 



- 9 - 



form of these functions significantly complicates the full evaluation of the backreaction, 
but exact solutions can be found in a particular limiting case. 

Taking t — > oo, corresponding to the universe at very late times, yields closed-form 
solutions. In this limit, the functions in (4.24) are heavily skewed towards a small krfe, and 
we can expand the normalisation functions and to zeroth order, 

Ak « -3, Bk « 0. (4.26) 

Solving (4.9-4.11) in this case gives the correction terms, for any m > 0, 

■D rjS-m (i + ^)(2 -I- m)(4 -Km) r(m)sin(7rm/2)' ^' ^ 



^2 m _^ m)(3 + m)(4 -|- m)(6 -|- m) T(m) sm{irm/2) 

+ <), (4.28) 



(3 + m)(6 + m 
(T) _ „ P (2 - m)2^ TT 



V 



6 + m 
For small m these are 



(1 + j„)(2 -I- m)(4 + m){6 + m) T{m) sin(7rm/2) 

(4.29) 



^m^^^ p(T) 8 (T) o(^)^«i7^(^) f4 30) 



In contrast to both the situation in radiation domination and that of scalar perturbations 
at any time, all three corrections are of the same order of magnitude. Using (3.9), the 
effective energy density of the backreaction is now 



^' {l + m){2 + m){4: + m){6 + m)r{m)sm{Trm/2) 8 



which as in radiation domination is positive definite. The normalised effective pressure 
(3.10) is 



n(^) - ge + llm + m^) TT ^ _8 (T) 



(4.33) 



yl—m 

(1 -|- m)(2 + m)(3 + m)(4 + m)(6 + m) T{m) sin(7rm/2) 9 
To linear order, the effective equation of state (3.11a) is 

(T) 8 17m 



The backreaction produced by gravitational waves acts to accelerate the universe with 
an equation of state close to that of a cosmological constant and intriguingly close to the 

(T) 

current observational bounds [48]. Furthermore, since m > 0, w^g is slightly more negative 

(T) 

than the zeroth-order approximation w^q = —8/9. The direct impact on the Raychaudhuri 
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equation (3.11b) is 



-3/3r?" 



2i-m + 9m + 



m 



vr 



(1 + m)(2 + m)(3 + m)(4 + m)(6 + m) r(m) sin(7rm/2) 



10 + 9m + rn? ^{T) 



6 + 2m 



'eflf 



3 18 ' 



(4.34) 



This is positive definite and larger than the effective energy density by approximately a 
third. As the effective equation of state suggested, the tensor backreaction deep in matter 
domination therefore acts as an effective dark energy. The evolution of the effective energy 
density and acceleration, 



^7^^) oc 



(T) 



OC 7/ OC a 



m/2 



(4.35) 



is markedly slow. The emergence of this effective dark energy is entirely natural - we 
have introduced no new physics and assumed only the existence of primordial gravitational 
waves with a red spectral index. 

The different behaviour of the backreaction in the epochs of radiation and matter 
domination arises purely from the different nature of the expansion at these times - and, 
in particular, how it modifies the dynamical backreaction. After matter /radiation equality 
the single mode h oc jQ{x) is converted into the two matter modes until, at very late times, 
the second matter mode is again negligible. Early in matter domination we therefore expect 
the backreaction to change significantly before settling down at late times in the limiting 
case found above. 





Figure 1: The approximation scheme for t = 1. Left: the normalisations and Bk and their 
approximations. Right: the integrand of TZd and its approximation. 
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The full integrals in matter domination are extremely complicated to solve. We employ 
an approximation scheme wherein we stitch together sequences of Taylor series for low 
values of x with an averaged behaviour at higher x. For example, in and we Taylor- 
expand about X G {1/2,1,3/2,2,5/2}, and employ the substitutions 



sin^(x) 



^, cos2(x) 



2sin(2x) 



vr 
2^ 



(4.36) 



for X > 3. This converts the intractable trigonometric functions into piecewise combinations 
of polynomials which are easily solved. Figure 1 illustrates this approximation scheme for 

(T) 

Ak and Bk and the integrand for TZj~, . 

We solve the resulting equations for the model with m = 1/160 and the results are 
shown in Figure 2. The left panel shows the evolution of the correction terms TZx>, Vt> 
and Qx> across matter domination. Shortly after equality, the curvature correction and 
kinematic backreaction decay rapidly, while the dynamic backreaction remains approx- 
imately constant. From a time of r/ ~ lOr/e onwards (corresponding to a redshift of 
z ~ Ze/1000 ~ 100) the decay slows and the corrections tend towards the late-time limiting 
cases (4.27-4.29). Also plotted in the figure for t < 1 are the values evaluated for radiation 
domination. The backreaction terms over matter domination evolve smoothly between the 
two limiting cases, as expected. The right panel of figure 2 shows the evolution of the effec- 
tive equation of state of the backreaction, which likewise tends smoothly from its radiation 



(T) 

value of iLigg ^ 1/3 to the present value of 



(T) 
eflf 



-8/9. 



10 - 





Figure 2: The effective backreaction fluid as a function of t = ij/rje- Left: the corrections TZ\~, , 
V\, and , scaled by rf / (3. Right: the effective equation of state. For late times the backreac- 
tion acts as a dark energy. 

Primordial gravitational waves, therefore, naturally give rise to a backreaction that 
acts in radiation domination as an additional radiative fluid but late in matter domina- 



- 12 - 



tion evolves to become a dark energy. The properties of this effective fluid are close to 
those required for a dark energy by observation. However, the effective energy density is 
unfortunately low. Deep in radiation and matter domination respectively it is 

^Tiv < ve) = ^ (^(^*^r> ^^2iv » ve) = 1^ (|^(^*^r- (4-37) 

If At = rAs, with r m 1/20 and As ~ 2.5 x 10~^, these are 

^'S'iv < Ve) ^ 10-'\hvr, ^'SPiv » Ve) ^ lO-'^hv)"^. (4.38) 

Phrased differently, the backreaction from gravitational waves produces an error in the 
Hubble rate today of 

^ « 10-'{hrio)"'/\ (4.39) 
-no 

Since m ^ 1 the normalisation (A;*??)™/^ from the pivot wavenumber and conformal time 
will be of order unity. While the gravitational waves in such a universe certainly produce 
an acceleration that acts as a dark energy, it is unsurprisingly orders of magnitude below 
that which has been inferred from observation. Even so, the discovery of an accelerating 
volume-averaged expansion in a vanilla (FLRW+perturbations) system remains significant. 



5. Discussion 

In this paper we have considered the backreaction that arises from primordial gravitational 
waves. While our motivation has primarily been those arising from inflation, the analysis is 
limited only by the requirement that the gravitational waves exhibit a red spectrum. Our 
first aim was to consider the tensor backreaction deep in radiation domination and compare 
it to that from the scalar modes. The two contributions are quantitatively similar: both the 
scalar and the tensor backreactions have effective energy densities which are approximately 
constant and effective equations of state ~ 1/3, and the average spatial curvature TZd 
is the dominant contribution. However, while the scalar effective energy density is of order 
QeS ~ 10^^ — 10^^, the tensor effective energy density is two orders of magnitude smaller. 
This suppression arises mainly from the small tensor/scalar ratio r. This implies that it 
is unlikely that the tensor backreaction has ever dominated over the scalars. Moreover, 
while in the scalar case the kinematical backreaction Qx> is negligible and the dynamical 
backreaction Pd an order of magnitude smaller than Tlx>, in the tensor case the kinematical 
backreaction is almost equal to the curvature correction while the dynamical backreaction 
is negligible. 

Our second aim was to examine the changes that occur as the universe passes from 
radiation to matter domination and evolves towards the current epoch. The impact on the 
curvature correction and the kinematical backreaction is dramatic, both of them decay- 
ing rapidly in the matter era. The dynamical backreaction, in contrast, remains almost 
constant and indeed grows slightly. At a conformal time approximately ten times that 
of equality, the corrections tend towards asymptotic power laws ~ 7^^(2+'^it)^ The curva- 

(T) 

ture correction TZj-, remains dominant, with the kinematical backreaction decaying only 
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slightly faster. It is the dynamical backrcaction that generates the difference between the 
radiation and matter cases: it tends towards a value approximately nine-tenths that of the 
curvature correction. The inhomogeneous corrections arising from tensor modes deep in 
matter domination, in direct contrast to those from scalar s, are therefore all of equivalent 
magnitude. This then allows small changes in the corrections to produces differences in 

(T) 

^eff (3-lla), with the ratios between them producing the intriguing nature of the late-time 

tensor backreaction. 

The effective equation of state of the tensor backreaction tends for small ht towards 

(T) (T) 

w^q' —8/9 + 17nT/54. It passes through w^f^ = —1/3 at a conformal time rj Wije, 
corresponding to a redshift z 100 and approaches its asymptotic values at late times. 
The dark-cncrgctic nature of the field thus emerges naturally throughout matter domi- 
nation although it acts to accelerate the volume expansion from very early times. The 
effective energy density of the tensor backreaction is approximately constant and evolves 
with conformal time as Q,^^^ ~ ry""^ ~ z"^/^. It is, however, extremely low and suppressed 
from the backreaction in radiation domination by a factor of order ut- 

In a pure-dust universe, which would correspond to the universe in the far future 
(assuming A = and no other dark energy fluid) , the backreaction reaches the asymptotic 
values. The tensor backreaction continues to slowly grow and exhibit an equation of state 
close to that of a viable dark energy. The effective fluid cannot be the observed dark 
energy for any sensible model since the amplitude is unsurprisingly low, but it nonetheless 
acts to accelerate the volume expansion. Given that we can find an effective dark energy 
from gravitational waves, it would be interesting to consider the more recent sources of 
gravitational radiation that emerge during the epoch of structure formation. It would be 
immediately interesting to examine the waves sourced by first-order scalar modes; in matter 
domination, first-order scalars act as a dark matter while first-order tensors act as a dark 
energy. How the second-order tensors act is an entirely open question, but their impact 
could be equivalent to the first-order tensors and they are of immediate significant interest. 

We conclude that dark energy, with observationally reasonable equations of state, can 
arise from volume- averaging in pure relativity and in realistic systems, introducing no new 
physics beyond a slightly red distribution of primordial gravitational waves. While the 
energy scale of this "dark energy" is extremely low - approximately fourteen orders of 
magnitude lower than the dark energy which is observed - the significance remains: a 
natural dark energy can emerge in a vanilla FLRW-|-perturbations system. 
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